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ABSTRACT 

Real geological formations may exhibit 
resistivity anisotropy in two ways: micro-anisotropy 
and macro-anisotropy. Micro-anisotropy is 
intrinsically anisotropic because of the 
micrvstructure of the formation. However, 
macro-onisotnypy is often due to electrical and 
electromagnetic well togging methods can achieve 
only limited resolution of the resistivity layering, for 
example, we often have to consider a collection of 
many thin layers as one composite layer, which is 
then macro-anisotropic. Macro-anisotropy is also 
found in cases of a fractured formation, in this paper; 
it is assumed that the resistivity ts the same In ait 
horizontal directions, but is different in the vertical 
direction, i.e. a transversely isotropic layered model. 

The determination of resistivity anisotropy is 
desirables as it may indicate the presence of 
otherwise unresolved thin layers and fractured 
formations. From a hydrogeological point of view, 
these may severely influence the hydraulic flow 
pattern in the ground. Thin clay layers in an 
otherwise sandy formation will lower the vertical 
hydraulic conductivity considerably and will deflect 
infiltration, and thin sand and gravel layers in an 
otherwise clayey formation may be serve as fast 
hydraulic conduction channels for polluted water. 

Three-component induction well logging may 
be the best method to determine the resistivity 
anisotropy. However, the tool has stilt not used In 
China. All data are typical lateral and dual 
induction logs. Neither lateral -well logging methods 
nor inductive well logging methods alone can 
resohe the anisotropy of the formation. However, a 
joint inversion of lateral and inductive data makes 



that anisotropy be taken into account and it can also 
resolve the coefficient of anisotropy, thus 
contributing to a mom detailed description of the 
formation resistivity. 

In this paper, an analysis of the importance of 
taking anisotropy into account in inverse modeling 
is presented, and it is shown how the combined use 
of lateral and inductive logs can resolve the 
coefficient of anisotropy of a formation. Through a 
synthetic 2D model, we show that inductive methods 
will only sensitive to the horizontal resistivity of a 
layer, while the thickness is undistorted. That is to 
say, we can determine the horizontal resistivity and 
formation thickness by inductive methods, but we 
can not determine the vertical resistivity by 
inductive methods, thus we cannot determine the 
coefficient of anisotropy. However, apparent 
resistivity of lateral methods can be approximate as 
the geometry mean resistivity of horizontal and 
vertical resistivities. And apparent thickness of 
lateral methods ts the multiplication of anisotropy 
coefficient and real thickness. Therefore we can not 
determine any parameters alone by lateral logs. 
However, a joint inversion of data from lateral and 
inductive logs may determine three parameters: the 
coefficient of anisotropy, horizontal resistivity and 
formation thickness. Synthetic data show that the 
Joint inversion method is feasible. 

INTRODUCTION 

Resistivity logs are the most important methods to 
determine values of residual and movable 
hydrocarbon saturation, However, The effects of 
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hydrocarbon saturation. However, The effects of 
anisotropy make that the logs can not response the 
full and true information of formation resistivities. 
For example, Induction loggings only response the 
horizontal resistivites in vertical well. To get 
horizontal and vertical resistivities of formation 
ones need to joint different logs, THagiwara et a). 
(1999) identified and quantified resistivity 
anisotropy by joint induction and laterolog logs. 
However, we have found laterolog logs are almost 
only , sensitive to horizontal resistivity in vertical 
well since the function of the focus electrodes. 
Niels B. Christensen (1999) presnted the ability of 
the joint inversion of geoelectrical and transient 
sounding data to resolve macro-amsotropic layers 
and tested the results using realistic earth models 
determinded from electrical logs. Therefore, his 
paper enlightens us to use lateral logs. 

This paper presents the possibility to resolve 
resistivity anisotropy by joint induction and lateral 
logs. At the following, we first give theory basis, 
then discuss results of numerical simulations, finally, 
give some discussion and conclusion. 

THEORY 

To lateral logging, an anisotropic layer is 

equivalent to an isotropic layer of thickness , 

equal to the product of the coefficient of anisotropy 

and the true thickness, and of resistivity Pu?, 

equal to the mean resistivity (Keller and 
Frischknecht 196$), so that 

h IAT =2.h (l) 

PlAY =^-A, ~ yfp^Ph ~ Pm (2) 

Herein, p v ,p h is vertical and horizontal resistivity 
respectively. The coefficient of anisotropy A is 



defined by ^ = yjp v /p h , h is the true 
thickness. 

From the formulae (1) and (2) we cannot 
determine not only the coefficient of anisotropy X , 
but also the true thickness h. 

To induction logging, it only responses horizontal 
resistivity of a layer, while the thickness is 
undistoned, so that 

hi» D ~ h (3) 

Pum = Ph (4) 

From the formulae (3) and (4) we can get the true 

thickness h and horizontal resistivity p h of a 

layer directly when neglecting the borehole, should 
and invasion effect However, we still cannot 
determine the coefficient of anisotropy. Therefore, 
we cannot resolve the full resistivity information by 
induction or lateral logging separately. 
However, we can first resolve the true thickness 

h and horizontal resistivity p h by induction log. 

Then, to substitute h and p h into (2), we can 

resolve the coefficient of anisotropy /I . Thus we 
can get the horizontal resistivity, the coefficient of 
anisotropy and the thickness of each layer by joint 
iversion of induction and lateral logs. 

NUMERICAL SIMULATION 

First, we assume a 2D three-layer formation 
model. The parameters in each layer are labeled in 
Fig la. In Fig la, we only change Rt-v. Rt-v is 
set as 10,JOO,1000,and 1O000. The response results 
are shown in Fig lb. It is 0.5m lateral log. Wc 
can see that lateral logs have a certain extent 
sensitive to resistivity anisotropy. In addition, in 
fig 2a, we only change Rt-h. Rt-h is set as 
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I0,100,l000,and 10000. The response results are 
shown in fig 2b. W can «c that lateral logs are 
more sensitive to the change of horizontal resistivity. 
From fig lb and fig 2b, we can make a conclusion 
that lateral log responses the change of Rt-b and 
Rt-v at the meantime. Therefore, lateral log can be 
used to resolve the coefficient of anisotropy X 
when jointing with induction log. 

To present the possibility, wc firstly assume a 2D 
three -layer formation model. The parameters in each 
layer are labeled in Fig 3a. Next, we can get the 
synthesized lateral and induction logs by applying 
the relevant forward modeling codes. The 
synthesized curves are shown in fig 3b and 3c. 
The induction log is the response of standard dual 
induction tooL Therein, Rild is apparent resistivity 
of deep induction log and Rilm is apparent 
resistivity of middle induction log. 

First, we can read the formation position is 
4 .2 5 -6. 25m and the true thickness h =2m from fig 
3 b of induction logs directly. 

Next, the horizontal resistivity p k can be 

resolved by inversion of induction logs. The 
layering of Rilm can be chosen as the initial value of 

p h . The initial value of 3 layers is chosen as 2.1. 

1.08, 2.1 respectively. After 3 iterations wc 
almost get the true value of Rt-h. The result is 
shown in fig 4a. Inversion result has a good 
precision comparing with synthesized curve. 
Through in vera ion wfi can eliminate the effects of 
mud and shoulder. 

Then, we substitute h and p h into inversion 

process of lateral logs. We can resolve the 
coefficient of anisotropy X only after the 
inversion. The results are shown in fig 4b. We 
have got Rt-v as 8.5 only after 3 incerauons. Then 
X is 2.9154. Its true value is 3.162, The quality 
comparison result is shown in fig 4c. 



Finally, we have got the horizontal resistivity, the 
coefficient of anisotropy and the thickness of each 
layer by joint inversion of induction and lateral logs. 

DISCUSSION AND CONCLUSION 

The research shows that lateral logs have 
information about horizontal and vertical 

resistivities. In addition, induction log can provide 
the true thickness of formation and horizontal 
resistivity alone. Therefore, we can identify and 
quantify resistivity anisotropy by joint lateral and 
induction log. 

It is essential to inverse induction and lateral logs 
respectively to remove the effect of mud and 
shoulder. 

If there are multi-parameters to be inverted we 
have to use multi logs such as Rild and Rilm. 
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Fig 2a Formation model with change of Rt-h 
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Fig la Formation model with change of Rt-v 
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Fig 1 b Lateral log with change of Rt-v 
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Fig 2b Lateral log with change of Rt-h 
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Fig 3a Formation model 
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Fig 3b Induction log 



Fig 4b Inversion results of vertical resistivity Rt-v 
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Fig 3c Lateral log 



Fig 4c Lateral logs of different inversion iteration 
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